Implicit regularization beyond one loop order: scalar field theories 
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Implicit regularization (IR) has been shown as an useful momentum space tool for perturbative 
calculations in dimension specific theories, such as chiral gauge, topological and supersymmetric 
quantum field theoretical models at one loop level. In this paper, we aim at generalizing system- 
atically IR to be applicable beyond one loop order. We use a scalar field theory as an example 
and pave the way for the extension to quantum field theories which are richer from the symmetry 
content viewpoint. Particularly, we show that a natural (minimal) renormalization scheme emerges, 
in which the infinities displayed in terms of integrals in one internal momentum are subtracted, 
whereas infrared and ultraviolet modes do not mix and therefore leave no room for ambiguities. 
A systematic cancelation of the infrared divergences at any loop order takes place between the 
ultraviolet finite and divergent parts of the amplitude for non-exceptional momenta leaving, as a 
byproduct, a renormalization group scale. 

PACS numbers: ll.10.Gh, 11.15.Bt, 12.38.Bx 



INTRODUCTION 



The problem of defining a consistent perturba- 
tive regularization of dimension specific quantum 
field theories is by itself of theoretical interest. Chi- 
ral, topological and supersymmetric gauge theories 
belong to this class of theories. The very concept 
of chirality is highly dimension dependent as there 
are no Weyl fermions in odd dimensions. The 75 
matrix which enters the chiral projection opera- 
tors has no straightforward generalization to ar- 
bitrary dimensions. Hence, care must be exer- 
cised in applying dimensional regularization meth- 
ods. Yet the standard model as a full theory is 
actually anomaly free, from the practical stand- 
point perturbative calculations are necessary to 
test the standard model and extensions against the 
experiments. The appearance of spurious anoma- 
lies, which spoils renormalizability by violating chi- 
ral symmetry, demands the addition of symmetry 
restoring counterterms, order by order, in pertur- 
bation theory. This renders practical calculations 
very tedious. Besides dimensional regularization 
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(DREG) explicitly breaks supersymmetry because 
the number of degrees of freedom of gauge bosons 
and gauginos does not match for d 7^ 4. There- 
fore, an invariant regularization scheme becomes 
important in this case. An alternative approach 
based on algebraic renormalization, which discards 
an invariant regularization scheme, is based on the 
construction of Slavnov- Taylor identities which en- 
codes all the essential symmetries of the model in 
consideration, gauge and supersymmetry included. 
Such constraints, in turn, deliver unambiguously 
the (non-invariant) counterterms to be added in 
order to heal the symmetry breakdowns. The main 
drawback of this approach is that it is much sim- 
pler to work within an invariant scheme which 
hitherto is missing at least to all orders in per- 
turbation theory. Moreover, it is not a priori obvi- 
ous that, in general, there are no supersymmetry 
anomalies [l[ , although the evidences quoted so far 
points to supersymmetry being a full symmetry of 
the quantum theory |2j. 

From the phenomenological viewpoint, it is ex- 
pected that the minimal supersymmetric standard 
model (MSSM) can be realized at weak scale in the 
experiments at LHC around 1 TeV. This makes 
quantum effects important. That is to say, besides 
the direct detection of susy particles, supersymme- 
try can also be probed via the virtual effects of ad- 
ditional particles such as the eletroweak precision 
observables within the MSSM Examples are 
the correction for the W boson mass, the effective 
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leptonic weak mixing angle, the anomalous mag- 
netic moment of the muon. A sound perturbative 
calculation from the technical standpoint may also 
shed some light on the supersymmetry breaking 
mechanism by constraining the number of param- 
eters of the underlying soft breaking model when 
contrasted with experiments . 

Dimensional reduction (DRED), a framework 
where only momenta are treated in D dimensions, 
whereas 7 matrices and gauge fields remain or- 
dinary 4— vectors, is used in most practical cases 
where the supersymmetry breaking of DREG mat- 
ters, even though it has been proved mathemati- 
cally inconsistent Nonetheless, recent works 
have demonstrated that its validity can be ex- 
tended beyond one loop order by using the quan- 
tum action principle, although it remains to be cir- 
cumvented some problems related to QCD mass 
factorization From the pragmatical point of 
view, DRED is largely used to unravel supersym- 
metric extensions of the standard model [5| . How- 
ever it is important to note that in DRED it is 
unclear whether or to what extent susy is actually 
preserved. It is only known that some susy rela- 
tions are satisfied at the one and two loop level but 
even at one loop the checks do not exhaust all the 
Green's functions that could be affected by a susy 
breaking. 

An important coordinate space regularization 
independent framework, which operates in the 
physical dimension of the model, is Differential 
Renormalization (DR) . DR basically consists of re- 
placing coordinate-space amplitudes that are too 
singular to have a Fourier tranform by deriva- 
tives of less singular ones as well as an integra- 
tion by parts prescription. It has been proved 
to be a symmetry preserving scheme @. This 
is achieved by both constraining the number of 
independent scales which result from each diver- 
gence and applying a set of operational rules in 
the introduction of basic functions (Constrained- 
DR), which assume the same values as the orig- 
inal divergent ones except for coincident points 
(short distance behaviour). Whilst at one loop 
order Constraincd-DR automatically produces ex- 
pressions with a single mass scale which fulfill the 
corresponding Ward-Slavnov- Taylor identities of 
the underlying model, at two loops a consistent 
extension of this program has not yet been con- 
structed 0- 

Implicit regularization (IR) is a momen- 
tum space regularization independent framework, 
which, like Differential Renormalization, does not 
rely on dimensional extension of space-time dimen- 



sion (8j-[l9|. The basic idea behind IR is to as- 
sume implicitly the presence of a regularization 
as part of divergent amplitudes only in order to 
separate their regularization dependent from its fi- 
nite part, by using a simple algebraic identity in 
the integrand of the amplitude. This operation is 
analogous to apply Taylor operators in the inte- 
grand of the amplitudes in the BPHZ formalism 
[2l| . as far as mathematical rigourousness is con- 
cerned, with the advantage of not modifying the 
structure of the integrand. The divergences are 
singled out in terms of internal momentum inte- 
grals which need not be evaluated. Such proce- 
dure generalizes to higher loop order. These so 
called loop integrals may be cast solely as a func- 
tion of an arbitrary scale which both parametrises 
the freedom of separating the divergent part of an 
amplitude and plays the role of renormalization 
group scale. Accordingly, a minimal subtraction 
renormalization scheme in which such loop inte- 
grals are subtracted in the definition of countert- 
erms naturally emerges. An essential by-product 
of isolating the divergencies in IR are surface terms 
expressed by differences between loop integrals of 
the same superficial degree of divergence. This is 
a crucial point of IR. Although in principle they 
are arbitrary numbers, such surface terms can be 
shown to be related to momentum routing invari- 
ance in a Feynman diagram, namely the possibil- 
ity of a shift in the momentum integration vari- 
able. The surface terms play an essential role in 
preserving symmetries in Feynman diagram calcu- 
lations: they translate into finite local countert- 
erms, whose value is in principle arbitrary but can 
be determined by the symmetry requirements of 
the underlying theory. 

For an account on the connection between gauge 
invariance and surface terms in IR see for 
the study of chiral and gravitational anomalies and 
CPT violation see and [l6[ in which we also 
showed that IR is an ideal arena to display consis- 
tently the physics of models which have finite but 
undetermined parameters, as discussed by Jackiw 
in (2^ |. App lications to supersymmetry can be 
found in [14| . where the j3- function of the Wess- 
Zumino model is calculated to three loop order and 
in [ijj, for the calculation of the anomalous mag- 
netic moment of the lepton in supergravity. A gen- 
eralization to higher loop order was firstly sketched 
in [l3|. The role of the surface terms, as unde- 
termined parameters in effective phenomenological 
models, have been considered in [l8j]-[l!|. Along 
the lines of IR, a complete systematisation of one- 
loop four-dimensional Feynman integrals appears 
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in 0. 

The purpose of this work is to clarify some issues 
for calculations beyond one loop level within IR 
which are basic to any field theoretical model. In 
doing so, we pave the way for a complete systemati- 
zation of IR as an automatic invariant framework 
order by order in perturbation theory, i.e. that 
preserves gauge invariance and supersymmetry, as 
we have shown at one loop order. 

Those issues are: 1) Can we still display the di- 
vergencies as basic loop integrals in one internal 
momentum in any loop order, overlapping diver- 
gencies included? 2) In addition, may such basic 
loop integrals be written as a function of an arbi- 
trary parameter A which plays the role of renor- 
malization group scale in a general case whilst the 
subtractions are dictated by the BPHZ forest for- 
mula? 3) In connection to this, can the derivatives 
of the basic divergent integrals with respect to A 
also be displayed by loop integrals (or constants) 
in the lines of Implicit Regularization? 4) How do 
the surface terms, which encode momentum rout- 
ing invariance in the loops, look like in arbitrary 
loop order? 5) In dimensional methods ultraviolet 
and infrared divergences become mixed. How does 
IR deal with this problem, in general, for one loop 
and higher loop order amplitudes? 

In order to answer these questions we use the 
simplest renormalizable scalar field theory namely 
</> 3 theory in 6 dimensions as a working example 
and show how we calculate renormalization group 
functions in the IR framework. We also outline the 
IR rules to work out amplitudes to arbitrary loop 
order. 



II. IMPLICIT REGULARIZATION AND 

HIDDEN PARAMETERS IN BASIC 
DIVERGENT INTEGRALS AT ONE LOOP 
LEVEL 

In this section we state the basic steps of IR in 
one loop calculations. We also construct an arbi- 
trary parametrization of the divergences expressed 
by loop integrals in order to make contact with 
other regularizations. Also we show how the in- 
frared cutoff in the propagators cancels out in a 
subtle interplay between the divergent and finite 
parts of the amplitude at one loop level, the di- 
vergent part being expressed as a basic divergent 
integral independent of the external momenta. We 
deal with higher loop in the next section. 

1. In order to give mathematical rigor to any 
algebraic manipulation performed in the am- 



plitude, we implicitly assume that a regular- 
ization has been applied. It can be main- 
tained implicit, the only requirement be- 
ing that the integrand of the amplitude nor 
the dimension of the space-time is modi- 
fied. After eventually performing Dirac ma- 
trix traces, group index contractions etc., 
we cast the momentum-space amplitude as 
a combination of basic integrals. Typical ba- 
sic integrals are: 

1 1 Ifj,i Ifj,u, ■ ■ ■ = 

f d 2w k 1, k^, k^k v , • • • 

J {2n) 2w [k + h) 2 -m 2 }... [{k + k N ) 2 - m 2 ] ' 

for a ./V-point function in a 2w integer di- 
mensional space. The internal momenta ki 
are related to the external momenta Pi by 
Pn = h-k N , pi = k l+ i-h L , i = 1, .. .,N-1 
if we choose J2j=i Pj = 0- 

2. In the basic integrals, the divergent part 
is subtracted as basic divergent integrals 
which are obtained by applying recursively 
the identity 

1 _ 1 

(p — k) 2 — m 2 (k 2 — m 2 ) 

p 2 ~2p-k 

(k 2 — m 2 ) [{p — k) 2 — m 2 ] 

until the divergent part is free from the ex- 
ternal momentum (p) dependence in the de- 
nominator. This will assure local countert- 
erms. The basic divergent integrals have the 
general form 

f A d 2w k g^ H 
J (2ir) 2w (k 2 - m 2 Y ' 

/d k k^ 1 k ll2 ■ ■ ■ k^ n 
{2ir) 2w (k 2 -m 2 ) a ' 

where the superscript A indicates that the in- 
tegral is regularized and w — a > for the first 
integral and 2w — 2a + n > for the second. 
Eventual even powers of internal momenta in 
the numerator are simplified by adding and 
subtracting a mass squared term. 

3. Express the basic divergent integrals for 
which the internal momenta carry Lorentz 
indices as a function of surface terms (i.e. in- 
tegrals of a total divergence). For example, 
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in four dimensions: 



where b 



d 4 k 



(2tt) 4 (fc 2 - to 2 ) 3 



d 4 fc d 



+ 



4 J (2tt) 4 dk» \{k 2 ~m 2 ) 2 
<- A d 4 k 1 



(2tt) 4 (fc 2 - to 2 ) 2 ' 



Such surface terms are regularization depen- 
dent (e.g. they vanish in dimensional regu- 
larization). The possibility of making shifts 
in the loop momenta means that the sur- 
face terms are null which reveals momentum 
routing invariance (MRI) . A constrained ver- 
sion of IR (CIR) assumes that such surface 
terms are canceled by local (MRI) restor- 
ing counterterms. In practice, this is auto- 
matically realized by setting them to zero 
from the start. CIR is an invariant regu- 
larization which preserves the Ward-Slavnov- 
Taylor identities of the amplitude, save when 
quantum symmetry breaking occurs. In this 
case they should be considered as finite ar- 
bitrary parameters to be fixed on physical 
grounds. That is because anomalies are 
linked to momentum routing dependence in 
Feynman diagram calculations. 

4. The basic divergent integrals which encode 
the ultraviolet behavior of the amplitude 
need not be evaluated. We adopt the fol- 
lowing notation: 



hog{m 2 ) = 



A >2 



d 2w k 



1 



(2n) 2w (fc 2 - m 2 ) 1 



(2) 



Iquad{m 2 ) 



d 2w k 



(2tt) 2w (fc 2 - to 2 ) 



— ,etc. (3) 



for logarithmically and quadratically basic 
divergent integrals, etc. . 

5. These objects can be subtracted as they 
stand in the definition of renormalization 
functions through, for instance, the defini- 
tion of local counterterms in the process of 
renormalization. A minimal, mass indepen- 
dent scheme is defined by substituting m 2 
with A 2 ^ using the regularization indepen- 
dent relation in 2w dimensional space-time. 

hogim 2 ) = I log {\ 2 ) + b 2w {-l) w In (^) , (4) 
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?4^r^y> and subtracting 
hog{^ 2 ), A playing the role of renormaliza- 
tion group scale in the renormalization group 
equation. For infrared safe massless models 
a systematic cancelation of ln{m 2 ) steming 
from (j4|) and from ultraviolet finite part will 
render the amplitude well defined as m 2 —> 0. 

6. The remaining ultraviolet finite integrals are 
evaluated as usual, using Feynman parame- 
ters or an extensive library of methods in the 
momentum space (23[. 

A systematic presentation of the finite part of 
one loop A-point Green's functions in four dimen- 
sions is given in [2(3]. 



As a matter of illustration consider the self en- 
ergy graph of the massless (f>\ theory as an example, 
which is power counting quadratically divergent, 
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1 



2 J k (fc 2 — m 2 ) [{p — fc) 2 — to 2 



A 

From now on we adopt the notation J, = 

/ d 2w k/{2ir) 2w in 2w dimensions. An infrared cut 
off is introduced in the propagators and the limit 
to 2 — * will be taken at the end of the calculation. 
In order to separate the divergences, the identity 
(p} is applied three times to yield 



1 



+ 



P 



^{P 2 ) = % 

4(p • fc) 2 
"(fc 2 -m 2 ) 4 ' (fc 

(p 2 - 2p ■ fc) 3 
" (fc 2 - m 2 ) 4 [(p - fc) 2 - m 2 



P 



£ 1 (fc 2 — m 2 ) 2 (fc 2 — m 2 ) 3 

4 



,2^4 



(6) 



Note that the first three integrals on the r.h.s. 
are divergent and regularization dependent. The 
first one is quadratically whereas the other two 
logarithmically divergent integrals. The remain- 
ing integrals are ultraviolet (UV) finite and can be 
evaluated using Feynman parameters. However, in 
a massless case they are still infrared divergent. In 
the sense of IR, in which we do not evaluate loop 
integrals, we will show that a cancellation of the 
infrared divergences coming from the UV divergent 
and finite parts will always take place before taking 
the limit m 2 — » 0. 

Let us start by discussing the regularization de- 
pendent terms. The basic one-loop regularization 



5 



dependent objects of a six dimensional theory are 
obtained by choosing w = 3 in eq. |2]) and ([3]). 

In equation ([6|) there appears, besides the 
quadratic divergence which vanishes for a mass- 
less theory if a adequate parametrization is used, a 
difference between logarithmic divergent integrals. 
Let us write — zS — — iSoo — i^F, (— i£oo) — i^F 
standing for the power counting (divergent) finite 
basic integrals. Then according to the rules of IR 
we have 



= 2pV' 



fc (fc 2 -m 2 ) 4 2 log[m ' 



2p^{T%-- 9llv I log {m 2 )), (7) 



with 



" A (fc 2 - m 2 ) 4 



7" / 2 \ 



1 



a 



6 y fc dk» \(k 2 -m 2 ) 3 



(8) 



whereas £ f can be easily evaluated and has a sim- 
ple form in the limit m 2 { 0. If we write 

= ag^ (9) 

a being an arbitrary constant, ^ reads 



1ap z 
b 



6 



p In 



-rllogim ) + 



Ah 



(-—) 

V m 2 / 



9 



(10) 



with b = i/[2(47r) 3 ]. Using relation Q for six di- 
mensions, Iiog(m 2 ) = /; og (A 2 ) -bin , A 2 ^ 0. 
In the equation above, it is clear that the divergent 
logarithms as m 2 J, cancel out as we indicated 
before. We show in the next sections that such 
infrared divergence cancellation mechanism takes 
place for arbitrary ./V-point functions at one loop 
order as well as for higher loop order (see also [14[ ) . 
The final answer is 



i£(p 2 ) = - 



9 9 

P 9 



hog(>? 



in which we have taken a — 0. It is evident that the 
momentum routing invariance requirement a — is 
not essential for such simple scalar theory. In fact 



the choice of the routing in this case contributes 
as a finite local counterterm which is fixed in the 
definition of the renormalization scheme, or cquiv- 
alcntly in our case in a redefinition of the arbitrary 
constant A. Within IR, the ultraviolet behaviour 
is displayed in terms of basic divergent integrals 
such as Ii og (X 2 ) above, etc. . Their generalization 
to higher loop order as their (regularization de- 
pendent) explicit evaluation may contaminate the 
underlying physics by breaking symmetries and in- 
troducing junk constants. The derivatives of the 
basic divergent integrals with respect to the scale 
A (which are important for the calculation of renor- 
malization group functions) are either finite or may 
be written in the form of basic divergent integrals 
being therefore regularization independent as well. 
The simple subtraction of Ii og {\ 2 ) in the defini- 
tion of renormalization constants defines a mini- 
mal mass independent scheme in IR. In section Hvl 
we calculate the renormalization group functions. 

Before proceeding to ./V-point functions and 
higher loop order, let us make an interlude in or- 
der to illustrate the role of the surface terms, which 
encode a specific momentum routing in the loops 
as we discussed in the introduction, on gauge sym- 
metry. The generalization to higher loop order is 
subject of a forthcoming contribution. 

The one-loop QED vacuum polarization tensor 
in 4D, as it is calculated in [1(| with arbitrary 
momentum in the internal lines, k\ and k 2 : 

n M „ = U(p 2 )(p fJi p^ - p 2 g llv ) 



4 ( a\9^ - ^(kf + k 2 )a 2 g^u 



-(fc"fc 1 



k 2 ^2 



k"k^)a z g {llv g af}} 



(fci + k 2 ) a (k 1 + k 2 ) ll OL 2 g v 



- 2^? fc i + k%k$)g fiv a 2 g a f3 



In the equation above, p = k\ 
momentum and 



(12) 

k 2 is the external 



nV 



4 r 

3 



WA 2 ) + 



(4tt) 2 



In 





i 




3(4tt) 2 J 



includes the basic divergent integral. We have cho- 
sen the massless limit just for the sake of simplicity. 
Now the momentum routing dependent terms are 
proportional to c^'s, namely 



o^ig^v 



a 2 g^ v 



9^ 



9/j.u 



k (k 2 — m 2 ) 2 



(fc 2 — m 2 ) 2 ' 
k^ ki> 

Q 2^3 



(fc 2 



(13) 



(14) 



G 



and 
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1 



k (fc 2 — m 2 ) 2 
kfj,ki/ kakp 

2 _ ™2\4' 



(15) 



These parameters are surface terms. It can be eas- 
ily shown that 



d 



k dk^ \ (fc 2 — to 2 ) 2 / ' 



d 



(fc 2 — m 2 ) 



(16) 



(17) 



and 



r A d 


ki/k a 


L dk? 


(k 2 — m 2 ) 3 



9{^9a/3}{ct3 -a 2 ). 

In the case of gauge symmetry, both abelian and 
nonabelian, a constrained version of IR in which 
such surface terms are set to vanish delivers gauge 
invariant amplitudes automatically [15| . which is 
illustrated in the simple example above . The 
cancellation of these surface terms is also a re- 
quirement to preserve Supersymmetry 17J. That 
is to say: they represent the symmetry restor- 
ing local counterterms. In other words, momen- 
tum routing invariance seems to be the crucial 
property in a Feynman diagram in order to pre- 
serve symmetries. In fact such surface terms eval- 
uate to zero should we employ DREG to explicitly 
evaluate them. This property somewhat reveals 
why DREG is manisfestly gauge invariant yet it 
breaks supersymmetry because invariance of the 
action with respect to supersymmetry transforma- 
tions only holds in general for specific values of the 
space-time dimension. [29| . 

A particular situation, however, is the ocur- 
rence of quantum symmetry breakings (anoma- 
lies) . Anomalies within perturbation theory may 
present some oddities such as preserving a cer- 
tain symmetry at the expense of adopting a spe- 
cial momentum routing in a Feynman diagram e.g. 
in the (Adler-Bardeen-Bell-Jackiw) AVV triangle 
anomaly [25]. In the case of chiral anomalies, IR 
has been shown to preserve the democracy be- 
tween the vector and axial sectors of the Ward 
identities which is a good 'acid test' for regular- 
izations. The arbitrary parameter represented by 
the surface term remains undetermined and floats 
between the axial and vector sectors of the Ward 



identities. That is to say, in the anomalous am- 
plitudes, there is no possibility of restoring, at the 
same time, the axial and the vectorial Ward iden- 
tities. The counterterm that will restore one sym- 
metry causes the violation of the other and there- 
fore it does not make sense to set the surface terms 
to zero. The answer is to be established by phys- 
ical constraints on such amplitude. This feature 
has also been illustrated in the description of two- 
dimensional gravitational anomalies [16j |. 

We end this section by showing that we can 
parametrize the basic divergent integrals should we 
wish to make contact with an explicit regulariza- 
tion without assigning a definite value to a regu- 
larization dependent parameter hidden in an UV 
divergent amplitude. 

As we have argued before, in IR we need only 
the (loop) integral representation of the divergence 
and its derivatives. The latter can always be cast 
as a loop integral as well. Alternatively, we could 
construct a general parametrization for such ba- 
sic divergent integrals, which we do here for the 
purpose of illustration and comparison with other 
regularization methods. It is however not neces- 
sary for calculational purposes. 

For instance the basic logarithmically divergent 
integral at one loop order satisfies 



hog{m ) = —2, 



dm 2 



(19) 



with b — i/[2(47r) 3 ], from which we may construct 
a general parametrisation by integrating the equa- 
tion above. Namely 



A 2 

hogim 1 ) = -bln( — 



(20) 



where is an arbitrary constant and A is a 
mass parameter introduced for dimensional rea- 
sons, which dictates the UV behavior of the in- 
tegral. It is very important to notice that a 
free arbitrary parameter (3 appears. The explicit 
parametrization of the differences between loga- 
rithmically divergent integrals (fl~4"|) and (fT5"|) which 
represents the surface terms will render an arbi- 
trary constant in the end as the ultraviolet diver- 
gent logarithms are cancelled. 
Now, because 



d 

-g^Iquad(m 2 ) = 2Ii og (m 2 ), 
a possible parametrization of I qU ad{fn 2 ) is 



Iquad{m 2 ) 



(4tt) 3 



m4 







(21) 



(22) 
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Let us now evaluate the quadratic divergence in 
([3]). Recalling that In DREG, for dimensional rea- 
sons, we have a factor fi 2e , where fi is a mass pa- 
rameter, e = 6 — d, d is the space-time dimension, 
we have 



fJ> 2e Iquad(m 2 ) 



(4tt) 3 
IE + 0(e 



2 1 , , /-47T/i 2 

2 - + 1 + In ' 



(23) 



which vanishes as m? — > and resembles the 
parametrization (|22[) . Now, in order to compare 
IR with DREG let us calculate ([7]). Because 

47T^ 2 



H^h og {m z ) = -b 



e 



and that T°„ Dfl£;G = 0, we obtain 

2e^oo 2£ 1 2 r / 2n 



+ 7 B + 0(e) 
(24) 



i + ]n(-*=£] +7 * + 0( e ) 
e \ / 



,(25) 



which should be compared with l[7|). 



III. ONE LOOP N-POINT FUNCTIONS 



We now turn ourselves to a general one-loop am- 
plitude for a general massless amplitude, using the 
procedure described above. 

In a theory defined in d dimensions, with d even, 
one will frequently deal with logarithmically ultra- 
violet divergent integrals of the type, 

/ _ f A d 2w k 



1 



(2tt) 2 <" (k 2 - m 2 )[(fc + pi) 2 - m 2 ] 
1 



[(k + p n ) 2 — m 2 



(26) 



where n = w — I = d/2 — 1. The following expan- 
sion can be performed, 



I 



with 



1 



k (k 2 — m 2 ) n [(k + p n ) 2 — to 2 ] 4^ 1 



h = 



Pi + 2 Pl -k 



k (k 2 - m 2 ) l +![(/c + Pl ) 2 - m 2 ] 
1 



[(k + Pn) 2 — vn 2 



(27) 



(28) 



The expansion above has been obtained by using 
the identity (fTJ) in all factors containing external 
momenta pi dependence from i = l to i = n— 1. 
Each time the identity is applied, we obtain a new 
integral Ii that is ultraviolet finite. This is why we 
do not use the index A in these integrals. Next we 
want to show that in the sense of IR the limit m 2 — > 
is well defined. Our proof is done in two parts: 
first we show that the integrals Ii are well defined 
. A few algebraic manipulations which will lead to 
eq. (|32p will make that clear. Then we show how 
to handle the m 2 dependence in the regularized 
integral. 

We begin the calculation of by using Feyn- 
man's parametrization 



1 



ai • • • a. 



(l + a-2)\ 
(«-l)! 



dX 



(l-zi a^-i)"- 1 



[(ai - a t )xi + (o 2 - a l )x 2 H h a ; 



l+a-l 



where 
dX 



1 rl—xi 

dxi / dx2 
• < Jo 



1—331 



dxi-i. 
(29) 

In the integral Ii we have a — i + 1 and I = w—i+1, 
so that 



Ii = — / dA(l - a;i — 

(p 2 - 2p 2 i x w - l - 2{pi ■ p i+ i)x w -i-i 
1 



k (k 2 + Q 2 r+^ 



2{pi -p n )xi) 

(30) 



with 



Q 2 = p 2 n x\{\ — xi) H hp 2 x„_ i+ i(l - x„_ i+ i) 

-2 J^(p; • p t )x„_;+ia; n _t + i - to 2 . (31) 

The final result is given by 
i (—l) w f 

dX(l -x\- x w -d 

2(pi -Pn)xi) 

(32) 



(4tt)™ 11 

(p 2 - 2p 2 x w -i - 2(pi -pi +1 )x w -i-\ 
1 



It is clear from the expression above that whenever 
the external momenta are such that p 2 ^ 0, the 



limit 



is well defined. 



We are now left with the regularization depen- 
dent integral, 



I* 



1 



k {k 2 - TO 2 )"[(fc + Pn f 



(33) 
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from which we can separate the appropriate ba- 
sic divergence ifogim 2 ) (typical of this dimension) 
from its finite part: 



l k (k 2 - m 2 ) w [{k+p n ) 2 ~ m 2 } ' 

(34) 

with lf og (m 2 ) given by eq. ([2]), and where we must 
remember that n = w — 1. Note that the cutoff 
mass m 2 is present in both finite and divergent 
parts. We will now show that, due to a scale rela- 
tion, the m 2 dependence of the divergent integral 
can be extracted and precisely cancels out the m 2 
dependence of the finite part. First, we make use 
of the identity 



1 



(k 2 -m 2 Y {k 2 -\ 2 Y 



(A 2 - m 2 ) 



E 



(35) 



(fc2 _ m 2)i(jfc2 _ A 2)«,-i+l> 

and obtain the d- dimensional scale relation 

lf og (m 2 ) = lf og (X 2 ) + b d {-l) d ' 2 \n (^j, (36) 

with b d = (47r t )d/2 wam ■ Again, like in the six di- 
mensional case, the relation above could be ob- 
tained from a particular regularization , but, as 
we have seen, the relation is regularization inde- 
pendent. On the other hand, the calculation of 
the finite part yields: 

Pn ~ 2 Pn ■ k 



k (k 2 - m 2 ) w [(k+p n ) 2 - m 2 ] 

= {-i) d/2 b d J\x(i- X yi 2 - 2 

In (^\- X lr m \ (37) 



(— m 2 ) 

We can now clearly see, by collecting the two parts 
together, that 

I A = 4 9 (A 2 ) - (-l) d/2 b d g-l) J\x 

(l-xf/ 2 - 2 \n ^- 2) X) y (38) 



where the limit m 2 — > has been taken. The whole 
integral, 



(39) 



becomes then m 2 independent. This result, in con- 
junction with equation (|32|) . will be important for 
the generalization of IR to higher loop order in 
order to subtract subdivergences. Also, once the 
subtraction scheme is appropriately chosen, the 
object I? oq (X 2 ) can be subtracted without having 
to be explicitly evaluated. The calculation above 
establishes a procedure for dealing with massless 
theories at one loop order in the context of IR. 
This procedure can be straightforwardly adapted 
to higher orders when non-overlapping divergen- 
cies occur, as shown in the ref. [14j |. where the f3- 
function of the massless Wess-Zumino model was 
calculated at three loop order. In that contribu- 
tion, it was shown that a n-loop scale relation can 
be used in order to appropriately handle the de- 
pendence on the infrared cut-off. 

We shall now generalize this procedure for the 
general case when overlapping divergencies occurs 
with support in the BPHZ-forest formula. 



IV. HIGHER ORDER CALCULATIONS 



Renormalization is a recursive program. In a n- 
loop order calculation the main goal is to identify 
the typical divergence of the n*' 1 order and the finite 
part of an amplitude once the (n — l) th has been 
renormalized. In [14| . we have calculated within 
IR the three-loop (3 function for the Wess-Zumino 
model in which there were not overlapping diver- 
gences to consider. 

In fact IR can be applied in the sense of defin- 
ing loop integrals as basic ultraviolet divergent ob- 
jects also when overlapping divergences occur. A 
regularization independent scale relation will also 
emerge in this case which serves to both cancel 
infrared divergences in a similar fashion as we pre- 
sented in the last sections and introduce a renor- 
malization group scale. New surface terms appear 
beyond one loop level and are expected to play an 
essential role in preserving gauge symmetry. More- 
over IR rules are compatible with BPHZ forest for- 
mula, which judiciously defines the set of subtrac- 
tions to remove the subdivergences. The rules to 
implement IR beyond one loop order can be sim- 
ply stated once we adopt the version of the forest 
formula which is analogous to the ordinary coun- 
terterm method in which the subtraction operators 
are translated into a local counterterm which sub- 
stitutes the subgraph. This is implemented via a 
recursion equation which involves disjoint renor- 
malization parts only (2l| . described in texbooks, 
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for instance [2q . 

1. Starting from one loop order, the 
(sub)divergences should be expressed in 
terms of basic divergent integrals, which 
can be written in terms of one internal 
momentum only and an arbitrary non- 
vanishing parameter A, the renormalization 
group scale of the method. For this purpose 
identity (P) is judiciously used. We assume 
that subdiagrams are proper (one particle 
irreducible) which is sufficient to discuss 
renormalization in general. 

2. The counterterms are defined order by or- 
der through a minimal subtraction process 
within implicit regularization which amounts 
to subtracting the (sub)divergences as basic 
divergent integrals as defined above 

3. The subtraction of subdivergences follows 
the forest formula written in a way which is 
equivalent of the subtractions of local coun- 
terterms only. That is to say given a dia- 
gram G the subtraction operator for a cer- 
tain renormalization part H has the effect of 
crushing H to a point and multiplying the 
counterterm to the resulting Feynman inte- 
gral. In this recursive process only the set of 
disjoint renormalization parts is needed. 

Mathematically speaking, a recursion relation 
for the subtraction of subdivergencies based on the 
forest formula for a Feynman diagram G in which 
only the disjoint renormalization parts are needed 
can be written as 

RgFg = G/v 11 ( — ^RhFh)) 

= E II i-t H RH)F G (40) 

where Rq is the operation to subtract only sub- 
divergencies, Fq^h) is the part of the amplitude 
which represents the (sub)graph G(H), ip is a set 
of disjoint renormalization parts of G, namely 



divergent loop integrals as a function of the arbi- 
trary scale A. 

For the massless case the basic divergent inte- 
grals are very simple. For instance the logarithmic 
basic divergence to n-loop order is given by 

An)d ( 2, f A d 2w k 1 
Hoc ( m ) 



hi 



/ "' "" ' J (2tt) 2w {k 2 -m 2 ) w 
„ , , (k 2 -m 2 Y 



A 2 



(41) 



and the corresponding scale relation is written as 

Flog ( ) = Hog ( X )+2^ a l ln 



i=0 



! "') i n > ( ™_ 1 (42) 



A 2 



with the a[ w ^ depending on the dimension 2w. 



Let us consider the two-loop contribution to the 
self-energy in 0| theory. The nested subdiagram 
in fig. ([1]) can be calculated by simply substituting 
the subdiagram by its finite part. The remaining 
integral, in one internal momentum, includes only 
the divergence of the order and the finite part. 




FIG. 1: The 2-loop diagram with a nested subdiver- 
gence for the (f>e self-energy 

The amplitude for the diagram depicted in fig. 
dJ) reads 



2 J k k\k-p) 2 J l l 2 {l-kf 

(43) 

The integral in the momentum I has been calcu- 
lated before (eq. ([ITj)): 



h 2 ( ( P 



7 /o9 (A 2 ) + 61n 



A 2 / 3 



b\. (44) 



tP = {H/H dG,H = proper, disjoint, d H > 0} , We take itg finite part and obtain . 



G/i/j representing the diagram obtained from G by 
crushing all H in ip to a point. The counterterm 
graph so obtained can be constructed with the loop 
integrals characteristic of IR order by order. Of 
course a further operation 1 — to is required to the 
overall divergence. The natural (minimal) renor- 
malization scheme in IR is to subtract the basic 



( 2 )^ _ _.9*b 



6 J k k 2 (k — p) 7 



(45) 



10 



where 

/(a) = 
Notice that 



k k 2 (k- P ) 2 n { X 2 



(46) 



- i£?V) = (1~ V)), (47) 

where — t corresponds to the subtraction operator 
which has the effect of crushing the nested sub- 
divergence to a point and multiplying the coun- 
terterm, which has been calculated in the previous 
order, to the resulting Feynman integral. 

Now we evaluate 1^ . We introduce an infrared 
cutoff m 2 as we have done before to rewrite 



r(2) _ 



1 



k (fc 2 — m 2 )[(k — p) 2 — rn 2 



■In 



fc — mf 
-A 2 



(48) 

which is ultraviolet quadratically divergent. The 
identity expressed by eq. Jl} is applied three times 
to yield 




4(p- fc) 2 



(fc 2 -rn 2 ) 3 (k 2 -m 2 ) A {k 2 ~m 2 Y 
(p 2 - 2p ■ fc) 3 



(fc 2 - to 2 ) 4 \{p - kf 



(49) 



We define the basic two-loop divergent integrals in 
6 dimensions, 



, (/ 2 -m 2 ) 3 



1 / (t 2 -m 2 ) 
In rr; 



and 



I (2) (m z ) = I - 

quad\ ) ^ (/2_ m 2)2 



(50) 



In - 



(I 2 -to 2 ) 
A 2 



(51) 

Clearly lj^ ai i(m 2 ) vanishes for massless theories in 
a regularization independent fashion, whilst the 
divergent basic integrals with Lorentz indices can 
be written in function of a surface term: 



k a kp 



(fc 2 - m 2 ) 4 



In 



(fc 2 — m 2 ) 
A 2 



6 



d 



k dk@ \ (fc 2 — m 2 ) 3 



k a ln /_ (fc 2 - m 2 ) 



A 2 



(52) 



The surface term will be set to zero for the sake of 
momentum routing invariance. Thus we have 

I {2) = -J&m 2 ) + lp 2 h og (m 2 )+I {2 \ (53) 
/( 2 ) being the finite part which is given by 



f(2) 



hi 



k (fc 2 — TO 2 ) 4 

(p 2 - 2p ■ fc) 3 



k (fc 2 — m 2 ) 4 [(p — fc) 2 — m 2 ] 



A 2 
In 



-A 2 



The finite integrals above can be easily evaluated 
using the identity, 



In a = lim - (a e — 1) 



(55) 



which after Feynman parametrization in the limit 
where m 2 — > is given by 



f(2) b 2 /, f™ 2 
1 = T8 P ln A 2 - 



61n| | -16 

TO Z 



31n 2 (-PlUlllnf-^ 
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(56) 



We see the infrared divergence, parametrized by 
rn 2 , appearing in the divergent and in the finite 
part. In order to observe the cancellation of these 
two contributions, we use the one-loop (??) and 
the two-loop scale relation, 



3, (m 2 



m 
A 2 " 



(57) 



so that eq. (|5"3"|) becomes 

^{-4 2 i(A 2 ) + ^ 09 (A 2 ) 



(2) 



El 

3 



I 



--In 2 
2 

Now equations (|35"]) and 
ig A bp 2 



(4) 

















18 I ? " ;; 
27 



together give 

10 
3~ 



4 2 i(A 2 )-^o 9 (A 2 ) 



-2 ln H-A^J-Y fol H-^J + l^^ !,i 



Let us now consider the two-loop overlapping 
contribution for the self-energy in 0g theory, rep- 
resented by fig. It is given by 



11 




FIG. 2: The 2-loop overlapping self-energy 



-i^\p 2 ) = i 9 -A 



,4 /-A 



1 



2 a, k^{k-if{p~kf{ P -if 



.(60) 



To solve the overlapping integral A, it will be easier 
to rewrite it as 

A = —A\ + A-2 + surface term, 

by using the identity [27l ]: 

d f(k-l) a \_lf P , (p-/) 2 



{k-lf {k-lf 



k 2 (p - kf 



(61) 



with 



D = k 2 l 2 {k-lf(p-kf(p-l) 2 . (62) 



We disregard the surface term on momentum rout- 
ing invariance grounds to write 



and 



A l = 



Ao 



1 



k , k 2 ( P ~k)H 2 (l-kf 



k , k\p-kfl 2 {p~lf 



(63) 



(64) 



In order to reduce to these two integrals, we have 
performed shifts (surface terms are not consid- 
ered). The integral Ai is actually a product of 
two independent one-loop integrals in the internal 
momenta I and k. It gives us 

P 2 lr2 ,^2^ 14 



m = yf og (x 2 ) Y bi log (y) 

+2b\n ! ^ //„,•: A" ; 



(65) 



For the A\ , after the integral in I is solved, we have 



A x = - 

+b In 
1 



i r i_ 

3 J k JF~k) 
k 2 



-Aho g {\ 2 ) 



3 | hogi* 
1 



A 2 J 3 

2 ^ °6 



A 



1 



(p - ky 



In 



HA (p-fc) 4 
A 2 



(66) 



Using an infrared cutoff, m 2 , the first term above 
is proportional to I qU ad(jn 2 ), plus a surface term. 
Therefore it does not contribute in the massless 
limit if we call for momentum routing invariance. 
Hence 



A 



1 



3 Jk \{p ~ k) 2 — to 2 ] 2 

b r A i 

3 Jk [(P ~ k f - m 2 
1 



In - 



(k 2 - m 2 ) 
A 2 



p 2 - 2p ■ k (p 2 - 2p ■ k) 2 



(k 2 -m 2 ) {k 2 -m 2 f (k 2 -m 2 f 

(p 2 - 2p ■ kf 
(k 2 — m 2 ) 3 [(fc — p) 2 — m 2 



(67) 



It is easy to check that the first term above is 1^ 
If we use identity ((T|) in I^ 2 \ we find 



p 2 — 2p ■ k 



k (fc 2 — m 2 f[(p — k) 2 — to 2 ] 



hi 



(fc 2 - TO 2 ) 



A 2 



(68) 



and 



(p 2 - 2p ■ kf 



hi 



k (fc 2 — m 2 f[(p — fc) 2 — to 2 j 

= / (2) -0- 2 )+P 2 0- 2 ) 
= I^+p 2 I%(m 2 ). 



(fc 2 - TO 2 ) 



A 2 



(69) 



The last term of eq. (|6"?| is finite and can be cal- 
culated with the use of eq. (|55|) . The final result 
of Ai is 



Ai = ^p 2 { WA 2 ) + bin (-pj + 7 -b\ , (70) 
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and we have, for the overlapping diagram, 



at the g 4 order is given by 



-^V)=^{-4 2 0fl (A 2 ) 
+ ^bI log (X 2 )-2b\n(-^\l log (X 2 ) 



b 2 \n 2 { - 



P 



2 x 34 
--b 2 
9 



(71) 



The result above refers to the complete diagram, 
which includes the subdivergences with non-local 
terms. The diagrams, according to the forest for- 
mula, for the counterterms to be added so as to 
cancel the subdivergences are given in fig. (|3|). 





(a) 



00 



FIG. 3: Counterterms to cancel the subdivergences of 
the overlapping diagram 

We have two equal contributions that furnish us 



-i^cLub = -2ig 2 I log (X 2 )(-i^(p 2 )) 
= ^{/f 0ff (A 2 ) + 61n(-^)l iog (A^ 



-bh og {\ 2 ) 



(72) 



After adding the above counterterm, we finally ob- 
tain 



^n 2 (- P - 



e(1-( w -( (j) )H^(p 2 )), 



(73) 



where tr a ) and t^) are the subtraction operators 
corresponding to the renormalization parts in fig. 

As in the one loop calculation, the limit m 2 — » 
was taken in the end and we made use of a scale re- 
lation given by eq. (|57p. The complete self-energy 



*£ (2 V) = ^U 2 09 (A 2 ) + 



+£"2(A a )-^&WA a ) 



3 log 
-ffi<V), 



(74) 



where the tilde refers to the finite part. We are 
now able to calculate the anomalous dimension (7- 
function). As usual, we define the renormalization 
constants, Z3 and Z gi so that 

4>o = Zl /2 cj) and g = Z g g, (75) 
and the Lagrangian is written as 



2 v ^ 3V 2 



(76) 



The expansion in terms of the coupling constant, 
g, is given by 



Z 3 = l + a l9 2 + a 2 g 4 + 0(g b ) (77) 



and 



Z g = l + Pl g 2 + p 2 g 4 + 0(g b ). (78) 
According to our calculations, we have 

(79) 



"1 = ~7:hog{>?)i 



+ y^o 3 (A 2 ) 



(80) 



and 



zl /2 Z g = \-ig 2 I log {\ 2 ), (81) 
which gives us 



3 

Pi = -j^o 9 (A )• 
The 7-function is defined as 

A dZ 3 Z 3 
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2Z a dX 2Z a 



(82) 



(83) 
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where 



V. CONCLUSIONS 



f-A 9/ -2A 2 9f 

f - X dx~ 2X d>?- 



(84) 



So, we have that 

7 = i {dig 2 + (2/3iai + a 2 - a^g 4 } + 0{g 6 ), 

with (3i the coefficient of g 3 in the /3-function, de- 
fined as 



Taking into account that 



(85) 



(86) 



^l> 2 )) = M^-WA 2 )] (ST) 



and 



we find 

and 

7 



^(^o S (A 2 ))- A. 



4(4tt) 3 



(88) 



(89) 



1 2 11 1 4 



12(47r) 3 * 108 4(47r)< 



g* + 0(g b ). (90) 



Notice that /3i and the g 2 contribution to the 7- 
function are universal results whereas the g 4 coef- 
ficient of the 7-function is renormalization scheme 
dependent. These results corresponds to the min- 
imal subtration scheme in implicit regularization. 



In this paper we have addressed important is- 
sues regarding the extension of implicit regulariza- 
tion beyond one loop order. Namely we showed 
that we can display the divergencies as basic loop 
integrals in one internal momentum in any loop or- 
der including the case of overlapping divergencies. 
This is the heart of IR which enables us to work 
in the physical dimension of the model. Moreover 
we have shown the general form of surface terms 
which appear as finite differences of divergent in- 
tegrals with the same superficial degree of diver- 
gence. They are important because they appear 
to relate momentum routing invariance to gauge 
and supersymmetry in a n-loop amplitude. We 
have explicitly verified this link in one loop order 
for both gauge and supersymmetry and to three 
loop order for supersymmetry. A general proof of 
the connection between momentum routing invari- 
ance in a general Feynman diagram and the vital 
symmetries to be respected by the corresponding 
Green's function is interesting to be constructed 
[28l | . Besides We have also shown that basic n-loop 
integrals may be written as a function of an arbi- 
trary parameter A which plays the role of renor- 
malization group scale; they are absorbed as they 
stand in the definition of the renormalization con- 
stants in the light of the BPHZ forest formula in 
a regularization independent fashion. That is be- 
cause the derivatives of the basic divergent inte- 
grals with respect to A can be displayed by loop 
integrals as well. Finally, contrarily to dimensional 
methods where ultraviolet and infrared divergences 
may become mixed and lead to ambiguities, IR 
clearly separates ultraviolet and infrared degrees of 
freedom through a kind of scale relation obbeyed 
by the basic divergent integrals. 
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